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A perturbative quantum master equation is derived for a system interacting with its environment, 
which is more general than the ones derived before. Our master equation takes into account the 
effect of the energy exchanges between the system and the environment and the conservation of 
energy in the finite total system. This master equation describes relaxation mechanisms in isolated 
nanoscopic quantum systems. In its most general form, this equation is non-Markovian and a 
Markovian version of it rules the long-time relaxation. We show that our equation reduces to the 
Redfield equation in the limit where the energy of the system does not affect the density of state of 
its environment. This master equation and the Redfield one are applied to a spin-environment model 
defined in terms of random matrices and compared with the solutions of the exact von Neumann 
equation. The comparison proves the necessity to allow energy exchange between the subsystem 
and the environment in order to correctly describe the relaxation in an isolated nanoscopic total 
system. 

PACS numbers: 05.50.+q, 03.65.Yz 



I. INTRODUCTION 



Studying the dynamics of a simple system interacting with its environment is a very important problem in physics. 
The theoretical description of this problem started a long time ago. 

In the context of classical mechanics several master equations, such as the Boltzmann equation, the Chapman- 
Kolmogorov master equation or the Fokker-Planck equation, were derived in order to describe the time evolution of 
the probability density of the system variables. 

In the context of quantum mechanics, wich interests us in this paper, the time evolution of a system interacting 
with its environment is described in terms of a reduced density matrix that is obtained by tracing out the degrees of 
freedom of the environment from the total (system plus environment) density matrix. In this way, the first quantum 
master equation was obtained by Pauli QyQ in 1928. This equation is called the Pauli equation and describes the 
evolution of the populations (i.e., the diagonal elements of the density matrix) when the system is weakly perturbed 
by an additional term in its Hamiltonian. The transition rates between populations are given by the Fermi golden 
rule. In 1957, Redfield p| derived the so-called Redfield equation in the context of NMR for a system such as a 
spin interacting with its environment. This equation has been widely used and applied to many systems where the 
dynamics of the environment is faster than the dynamics of the system. This equation is Markovian and has the defect 
of breaking the positivity on short time scales of the order of the environment correlation time for initial conditions 
near the border of the space of physically admissible density matrices. Many similar master equations for a system 
interacting with an environment have been derived since then starting from the von Neumann equation and making 
several assumptions (weak coupling limit, Markovianity, separation of time scale between system and environment) 
IB IE HI In 1976 Lindblad [l(j derived the most general quantum master equation which is Markovian and 
which preserves positivity. The Redfield equation has a Lindblad form in the case of 5-correlated environments. More 
recently, a non-Markovian Redfield equation has been obtained that preserves positivity and reduces to the Redfield 
equation in the Markovian limit It has also been shown 0, that the Markovian Redfield equation can 

preserve positivity if one applies a slippage of initial conditions that takes into account the non-Markovian effects on 
the early dynamics. Similar considerations have been proposed for different master equations |l3l Il4l Il5| . As far as 
one considers the weak-coupling regime, all the master equations derived till now in the literature at second order of 
perturbation theory can be deduced from the non-Markovian Redfield equation. 

The problem is that the non-Markovian Redfield equation as well as the other aforementioned master equations 
existing in the literature are based on the fundamental assumption that the environment does not feel the effect of 
the system. This assumption seems realistic for macroscopic environments but not in the case of nanoscopic isolated 
total systems in which the density of states of the environment can vary on an energy scale of the order of the system 
energy scale. Because nanoscopic physics is experimentally progressing very fast, we expect that such effects will 
become important and measurable in future applications. Alre ady, quantum dissipation is being envisaged on the 
nanoscale for applications such as spin d yna mics in quantum dots [l6j or isomerizations in atomic or molecular clusters 
in microcanonical statistical ensembles [lTj. Another possible application is the intramolecular energy relaxation in 
polyatomic molecules |lgj . 
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The aim of the present paper is to systematically derive from the von Neumann equation a master equation which 
takes into account the fact that the energy of the total system (system plus environment) is finite and constant and, 
therefore, that the energy distribution of the environment is affected by energy exchanges with the system. The 
aforementioned equations can be derived from our master equation, which thus appears to be very general. 

The plan of the paper is the following. In Sec. [H] we systematically derive our master equation and the non- 
Markovian Redfield equation from the von Neumann equation, by pcrfoming a second-order perturbative expansion 
in the coupling parameter (under the assumption of weak coupling) for general environments. Thereafter, we consider 
the Markovian limit in both cases. We also show how, in this limit and neglecting the coupling between the populations 
and the quantum coherences, our master equation reduces to a simple equation of Pauli type for the total system, 
taking into account the modifications of the energy distribution of the environment due to the energy exchanges 
with the system. Finally, we compare our master equation to the Redfield equation and discuss how the Redfield 
equation can be seen as a particular case of our master equation. In Sec. IIIII we apply our master equation and 
the Redfield equation to the case where the system is a two- level system interacting with a general environment. 
In Sec. IIVI we apply the master equations to the case where the system is a two-level system interacting with a 
complex environment (such as a classically chaotic or many-body environment) that is modeled by random matrices 
from a Gaussian orthogonal ensemble, which we call Gaussian orthogonal random matrices (GORM). In Sec. [V] we 
compare the solutions of the non-Markovian and Markovian master equations to the exact solutions of the complete 
von Neumann equation in the case of our spin-GORM model. Conclusions are finally drawn in Sec. IVII 



II. DERIVATION OF THE FUNDAMENTAL EQUATIONS 

The Hamiltonian of the total systems that we consider here is made of the sum of the system Hamiltonian Hs an d 
the environment Hamiltonian Hb plus a coupling term that has the form of the product of a system operator S and 
a environment operator B. The generalization to a coupling term of the form X^2 t SiBi is easy. The amplitude of 
the coupling term is determined by the coupling parameter A: 

Htot = H a + XV 

= H S + H B + XSB . (1) 

The eigenstates of H$, respectively of Hb, will be denoted by \s), respectively \b). The eigenvalues of Hs, respectively 
of Hb, will be denoted by E s , respectively Eb- Finally, the eigenstates of Htot will be denoted by \a) and its eigenvalues 
byE a . 

The evolution of the total density matrix is described by the von Neumann equation: 

f){t) = -i[Htot,p{t)]=Ctotp{t) , (2) 

where £tot is the so-called quantum Liouvillian or von Neumann operator of the total system. The interaction 
representations of the operators are given by 

V(t) = e iHoty e -iH t j 

B(t) = e i6Bt Be- iABt , 

S(t) = e tAst Se~ lAst . (3) 

In the interaction representation, the von Neumann equation becomes 

pi{t) = -i[\V(t),p x (t)] = £j(i)M*) , (4) 

with the interaction Liouvillian £/(t) = e~ Cot Cie Cat where e Cot A = e~ lHot Ae lHot , ^ e f ree LiouviHian £q = Cs+Cb = 
— i[Hs, •] — i[Hs, ■], and A is an arbitrary operator. The perturbative expression of the von Neumann equation in the 
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interaction representation is given to order A 2 by 



p/(t) = p(0) + / tfti£/(ii)/5(0) 



i) 



/ dt! / 1 dt 2 jC I (t 1 )jC I (t 2 )p(0) + 0(X 3 ) (5) 



n JO 

t 



p(0)+ / dTe- CoT de CoT p(G) 
Jo 



[ dT [ dre- CoT £/e £oT /:/e- CoT e £oT p(0) + 0(A 3 ) , 
Jo Jo 



(6) 



if we set T = t± and t = t\ — t 2 . Equation © is the starting point of all the derivations of a master equation in the 
weak coupling limit for a total system made of a system and its environment in mutual interaction. 

A. Our quantum master equation 

We now derive our master equation which is the central result of this paper. 

The main idea is to describe the time evolution in terms of quantities which are distributed over the energy of the 
environment. We thus define the following quantities in terms of which we intend to describe the properties of the 
system: 

P 3S , (e; t) =Tr p(t)W)(s\ 6(e- H B ) . (7) 

The diagonal element P ss (e;t) is the probability density to find the system in the state s while the environment has 
the energy e. The off-diagonal element P ss > (e; t) characterizes the density of the quantum coherence between the 
states s and s' , density which is distributed over the energy e of the environment. 
The matrix composed of the elements P ss r (e; t) is Hcrmitian 

P ss ,(e;t) = P;, s (e;t). (8) 

Moreover, the normalization Trp(t) = 1 of the total density matrix implies that 

]T f de P aa (e;t) = 1 . (9) 

s •* 

In order to obtain a closed description in terms of the quantities Q , we suppose that the total density matrix can 
be described at all times by a density matrix of the following form: 

^)=£i*><^^#^ ( 10 ) 

n(H B ) 

where we have defined the energy density 

n(e) = Tr B S{e - H B ), (11) 

which is supposed to be smoothened on the energy scale of the mean level spacing. The assumption I|1Q(I has the 
effect of neglecting the contributions from the environment coherences to the system dynamics (albeit the system 
coherences are kept in the description). We remark that the form (|10|l is not supposed to strictly hold at all times 
but is an assumption in order to obtain a closed set of equations for the quantities P ss > (e; t). 

In order to better understand the meaning of the above definitions, we notice that the reduced density matrix of 
the system takes the form 

p s {t) = Tr B p(t) = J deTr B 5(e - H B )p{t) = J de^ \s)P ss >(e;t)(s'\ , (12) 
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which can be represented in the basis of the eigenstates of the system Hamiltonian as 



hit) = / de 



( P n {e;t) P 12 {e;t) ... Pm s {e;t) \ 
P 2 i(e:t) P 22 (e;t) ... P 2Ns (e;t) 



\PN s l(e;t) P Ns2 (e;t) ... P NsNs (e;t) J 



(13) 



The goal of the precedent choice for the form of p(t) is thus to obtain a description in which the state s of the system 
is correlated with the energy e of the environment. In other words, the density matrix ps of the system is decomposed 
as a distribution over the energy e of the environment. 

We now proceed to the derivation of the equations of motion for our quantities P ss >(e; t) in the weak-coupling limit. 
We start from the perturbative expansion © of the total density matrix in the interaction representation ©. We 
first define the interaction representation of our quantities Q): 

P Issl {e-t)^e^- E ^ P ss ,{e-t) . (14) 

We now have that 

P Iss ,(e;t) = Trp I {t)\s'){s\ S(e - H B ) . (15) 
Inserting the perturbative expansion ||SJ), we get 

P Isa ,(e;t) = Tri"p(0)+ / dT TiXe- c ° T Cie c ° T p{0) 

Jo 

+ f dT f dr TrXe- CoT '£ ie CoT £ ie - CoT e CoT ' p(Q) + 0(A 3 ) , 
Jo Jo 



(16) 



where X — \s')(s\ S(e — Hb). Differentiating with respect to time, we obtain the equation 

P Iss >{e;t) = TrXe- Cot £ie Cot p{0) 



+ [ dTTrle~ £ot £/e £oT /:/e- £or e £ot p(0) + 0(A 3 ) , (17) 
Jo 

where the initial density matrix takes the assumed form (|10|) with t = 0. 
The first term is thus explicitly given by 

TrXe- Cot £ie Cot p(0) = -i\Tr\s'){s\ 6{e - H B )e lkot [V, e~ l " ot p(0)e l " ot } e - l6ot 
= -age i ^- E '-) t (s|,S|s)P^( e ;0)n(e)(5) £ 

S 

+*A^ e -' i ^'-^ t (s|5| S ') J P ;5S -(e;0)n(e)< J B) e , (18) 

8 

with the environment coupling operator B averaged over the microcanonical state of the environment 

n(e) 

We now assume that this average vanishes, (B) e = 0. Otherwise, the nonvanishing average is absorbed in the 
system Hamiltonian by the following substitutions: 

H s -» H s + \(B) e S, (20) 
V -» V-(B)J, (21) 
H B ^ Hb , (22) 

leaving unchanged the total Hamiltonian. Thanks to this simplification, the first-order term of the perturbative 
expansion vanishes 

TrXe- Cot £ ie Cot p{0) = . (23) 
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As a consequence, the time evolution of the total density matrix is given by the uncoupled Hamiltonian Ho = Hs + H B 
up to correction of the order of A 2 : 



(24) 



According to our closure assumption that the total density matrix keeps the form l|lUfl during its time evolution, we 
have that 



_,. P S ff{H B \t) 



e^V(O) = P(t) + 0(A 2 ) = J2 \s){s'\ tlE^fl + (A 2 ) , 



n(H B ) 



(25) 



which we can substitute in Eq. I|17fl to get 

Pi ss >(e;t) = f dTTrA > e- £ot £/e £or /:/e^ £oT e £ot j 5(0) + 0(A 3 
Jo 

= -A 2 e^- B =')*XJ^^Tr{| S ')(s| S(e-H B ) 



SB, 



0(A 3 ). 



n(H B ) 

Going back to the original representation with Eq. I|14f) and expanding the two commutators, we obtain: 



(26) 



P ss /(e;t) = —i(E s — E S >)P SS > (e; t) — A 2 ^ J dr { 



(s\SS{-T)\s}(s'\s'}Tr B 6{e - H b )BB{-t) 



Pss>(H B ;t) 
n{H B ) 



-{s\§\s){s'\§(-T)\s')Tr B 6{e - ff.) /^'^ ;t) B(-r) 

n(H B ) 

-<.s|^(-r)| S -)( S -^| S 0Tr B 5( e -^ s )B(-r)^|^B 

n(H B ) 

+ ( S | S -)( S -'|5(-r)5| S 0Tr B 5( e -i? B )^fT^^(- T )^i +°( A3 ) ■ 

n(H B ) J 

In order to evaluate the four last terms, we notice that, for a quasicontinuous energy spectrum, we can write 

Tr B S(e-H B )BS(e'-H B )B = ]T S(e - E b )S(e' - E b ,)\(b\B\b')\ 2 

b,b' 



(27) 



= n(e)n(e')F(e,e') 



where the function F(e, e') stands for 



F(e,e') = "|(e|% 



l\\2v 



(28) 



(29) 



where |e) denotes the eigenstate \b) of the environment Hamiltonian H B corresponding to the energy eigenvalue 
E b = e. Equation Ij29(l supposes some smoothening of the squares |(e|i?|e')| 2 of the matrix elements of B over a dense 
spectrum of eigenvalues around the energies e and e'. The function l)29|) has the symmetry 



F(e,e')=F(e',e) 



(30) 



With the definition l|29|l and the identity 



de'6(e - H B ) = I , 



(31) 
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we can now write that 



Tr B 5(e - H B )BB(-r) Pii '^ B '' *> = P m ( e ;t) [ de'n{e')F^e')e+^-^\ 

n{H B ) J 

Tr B 5{e - H B )B ^fffgj *} g(_ T ) = n ( e ) / d e ^ 5 ,( e ' ;t )F(e, e'^^"^, 
n{H B ) J 

Tr B 5{e - H b )B(~t) Psb '&b; Q £ = n{e) f de >p s _ l{e >. t)F ^ e ') e -^-e')* j 



(32) 



(33) 



(34) 



(35) 



Tr B( 5(e - g B ) ^ *) B(-r)B = P^(e;t) / de'n(e')F(e, ^e^-^, 
n{H B ) J 

Accordingly, our quantum master equation finally takes the closed form 

P ss >(e;t) = -i(E s -E a ,)P ss ,(e;t)-\ 2 J2 J de'F(e,e') J dr { 

s,s' 

(s\S\s') (s'\S\s) P Sa ,(e;t) n(e') e +i(^'+E s -E 3 ,) T 
-(s\S\s) (s'\S\s') P- s - sl {t'-t) n(e) e +i(—'+^-E s ,)r 
-{s\S\s) (s'\S\s') Pss'(e';t) n(e) e -<^+B.-E,)r 

+{s'\S\s) (s\S\s') P sS ,{e;t) n(e') e -<(«-«'+^-B,) T j + (36) 

Equation (|36|l determines the time evolution of the distribution functions P ss i (e; t) describing the populations 
and quantum coherences of a system influencing its environment and is the central result of this paper. It is a 
non-Markovian equation because of the presence of the time integral in the right-hand side. 

In Eq. H36fl . the function n(e)F(e, e') determines the properties of the coupling to the environment and, in particular, 
the time scale of the environment. If this time scale is supposed to be shorter then the system time scales { B 2 -e } ' 
we can perform a Markovian approximation in Eq. I|36|) . Such an approximation is justified for a process evolving on 
time scales larger than the environment time scale. The Markovian approximation consists in taking the limit where 
the upper bound of the time integral goes to infinity and using the following relations: 



f 

Jo 



i 

d T e ±iujT = ±iV-+irS(cj) , (37) 

UJ 



where V denotes the principal part. 
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We finally obtain the Markovian version of our quantum master equation (|36|) as 

P ss '(e;t) = -i(E s -E s ,)P ss ,(e;t) 

s,s' 

x [(s\S\s') (s'\S\s) P Ss >(e;t) - (s\S\s') (s'\S\s') P„-(e;t)" 



-(s\S\s) (s'\S\s>) n(e) 
x / de' F(e,e') P ss ,(e';t) 



-V- 



e-e' + E s - E s 



e-e' + E s , - Eg, 
- 7T A 2 ^ { n(e + E, - E- s , ) F{e, e + E- s - Eg, ) 

s,s' 

x [(s\S\s') (s'\S\s) Pg s ,(e;t) + (s\S\s') (s'\S\s') P sS (e;t) 
~{s\S\s) (s'\S\s') n(e) 

x [F{e, e + E s ,- Eg,) P ss > (e + E s , - E» ; t) + F(e, e + E s - Eg) P ss > (e + E s - E s ; t)] } 
+<3(A 3 ) . (38) 

We notice that the use of this Markovian equation may require a slippage of initial conditions as shown in Refs. 
[HD3. In Eq. (EEJ), the last terms in 7rA 2 typically describe the relaxation to a stationary solution. The terms in 
i\ 2 modify the frequencies of oscillations and include the so-called Lamb shifts of the zeroth-order energy eigenvalues. 
Indeed, if we consider only the evolution of the off-diagonal matrix element P ss > (e; t) by neglecting its coupling to all 
the other matrix elements, we obtain the equation 



P ss >(e;t) ~ {-i 
with the energies modified by the Lamb shifts 



E s (e)-E s ,(e) 



-L, 



<(e)} Pss>(e;t) 



E s (e) = E s + A 2 Y^\{s\S\s) | 2 J de 1 n(e') F(e, e') V ^^^ 
and the damping rates 



-Eg 



0(A 3 ) 



(39) 



(40) 



ttX 2 Y, [\(s\S\s)\ 2 n(e + E s -E s )F{e,e + E s ~E- s ) 

8&S) 

+ |( S '|^|s)| 2 n(e + E s ,-Eg) F(e,e + E s , - Eg) 
+n A 2 {(s\S\s) - ( S '|^| S ')) 2 n(e) F(e, e) + 0(X 3 ) , 



(41) 



in agreement with the results of Ref. . We notice that the complete equations for the off-diagonal matrix elements 
couple in general different energies because of the integrals over the environment energy e'. 

The evolution equations for the populations of the states \s) of the system can be obtained by neglecting the 
contributions from the quantum coherences, i.e., by neglecting the terms involving off-diagonal elements of P ss i(e;t). 
This is justified in the weak-coupling limit as long as the coherences vanish or are negligible in the initial conditions, 
i.e., P ss '(e; 0) = for s ^ s' . Accordingly, we obtain the following evolution equations for the populations: 

P ss (e;t) ~ 2ttX 2 Y / \( s \SW)\ 2 F(e,e + E S ~E S> ) 

s' 

x [n(e) P sV (e + E a - E s , ;t) - n(e + E s - E s ,) P ss (e; t)] . 

(42) 

This equation is a kind of Pauli equation established with the Fermi golden rule and the conversation of energy in 
the transitions. Indeed, if a transition happens from a state in which the energy of the system is E s and the one 
of the environment e to a state in which the system has energy E s i , the final energy of the environment should be 
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e' = e + E s — E s i, which is well expressed by Eq. I|42|) . Nevertheless, Eq. I|42(l rules the populations of the states 
s of the system with the extra information given by the distribution over the environment energy e, which is not a 
feature of the standard Pauli equation and which turns out to be of importance to understand the relaxation inside 
a nanoscopic isolated system. 

Our Markovian master equation l|38() is more general than an equation for the populations because it also describes 
the time evolution of the distributions of the quantum coherences over the energy of the environment. 



B. Comparaison with the Redfleld master equation 

We now discuss the conceptual differences between our quantum master equation and another one known as the 
Rcdhcld master equation. This equation is well known in the context of nuclear magnetic resonance (NMR) where it 
describes the time evolution of nuclear spins interacting with their environment. 

The Redfield master equation describes the time evolution of the system density matrix obtained tracing out from 
the total density matrix the degrees of freedom of the environment 

ps(t) = Tr B p(t) . (43) 

The Redfield equation is derived by using the closure approximation that the total density matrix keeps the form 

p(t) = Ps {t) ® p B (44) 

during the whole time evolution, where pB does not depend on time. The Redfield master equation is derived in the 
weak-coupling limit by a method similar to the one of the previous section to get 

k(t) = -i[H s ,ps(t)} 

—\ 2 S { dTa{T)S(-T)p s (t) 



+X 2 Sp s {t) / dTa*(T)S(-r) 
Jo 

+A 2 [ dra(T)S(~T)p s (t)S 
Jo 

~\ 2 p s (t) [ dTa*(T)S(-r)S + 0(\ 3 ) , (45) 



with the correlation function of the environment operators 

a(t) = (B(t)B(O)) = Tr B p B B{t)B(0) . (46) 

Equation (|45|l is a non-Markovian Redfield equation. The non-Markovianity comes from the fact that the integrals over 
expressions containing the correlation function depend on time. The density matrix of the system can be represented 
in the basis of the system eigenstates as 

Vss ,(t) = (s\p s (t)W) . (47) 
In this representation, the non-Markovian Redfield equation has the following form: 



Vss'(t) 



i(E s - E s .) r) ss , (t) - A 2 J de'n(e')F(e, e') jf dr { 



S\S\S') (S'\S\S) 7fey(t) 



,+i(e-e'+E s -E s >)T 



(S\S\S) (S'\S\ S ') Vm (t) e +i(e-e+E 3 -E a )r 
(S\S\S) (S'\S\S') Vm (t) e -<( e - e '+^-B,,)r 



+ {S'\S\S) (S\S\ S ') Vss ,(t) e -i(e-e'+E s ,-E s )r j 

+0(A 3 ) . (48) 
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If the environment is large enough, the correlation function in Eq. I|45[l goes to zero after a certain time. This time, 
called the environment correlation time r corr , determines the time scale of the environment dynamics. If we perform 
the Markovian approximation that consists of putting the upper bound of the time integral in the non-Markovian 
Rcdficld equation to infinity, one gets the standard Redfield equation. We notice that, in doing so, the time evolution 
may be spoiled on a time scale of order r corr unless some use is made of some slipped initial conditions [Til Il2| . 
Performing this Markovian approximation, one gets the standard (Markovian) Redfield equation given by Eq. I(45|l 
with replaced by J °°. As shown in Refs. 0, 0], the use of this Redfield Markovian equation needs to be 
supplemented by a slippage of initial conditions. 

In order to compare the Redfield equation with our master equation derived in the previous section, we consider 
the case where the environment is initially in the microcanonical state: 



PB 



6(e - H B ) _ 6(e - H B ) 



Tr B £(e - H B ) 



n(e) 



(49) 



Having chosen the microcanonical density matrix 149(1 for the environment, the correlation function (|46j) takes the 
form 



1 



a(j, e) 



n(e) 
1 



Tr B 6(e- H b )B(t)B 

de' TT B 5(e - H B )B(T)8{e' - H B )B 



n{e) 

de'n(e') F(e,e') e i(£ " £ ' )r 



(50) 



In the basis of the system eigenstates, the Redfield equation takes the form 
Vss'(t) = ~i(E s - E 8 ,) r) ss >(t) 



- i\ 



de' n(e') F(e,e') V- 



1 



^ { / < - ( ' - £\, - E 

\{s\S\s') {s'\S\s) r, Ss /(t) - (s\S\s') (s'\S\s'} r, sS (t) 

1 



Ee> — Ef, 



-V- 



E M - E K 



+ J de' n{e')F(e,e') 
x (s\S\s) (s'\S\ S ') Vss ,(t) } 
- vr A 2 { n(e + Eg - Eg,) F(e, e + E- s - E w ) 

x [(s\S\s') (s'\S\s) Vss ,(t) + (s\S\s') (s'\S\s') r,, s (tj 

- [n(e + E? - Eg,) F(e, e + Eg, - E s ,) + n(e + E s - E s ) F{e, e + E s - E.)] 
x (s\S\s) (s'\S\s') r, S g,(t) } 
+0(A 3 ) . 

The off-diagonal elements of the system density matrix individually obey the equations 



lis 



(t) ~ { -i [E s (e) - Eg, (e)] - T ss , (e) } Vss , (t) , 



(51) 



(52) 



with the same Lamb shifts (|4Uf) and damping rates H41|l as in our master equation and as expected from Ref. 0. 
There is no difference between our quantum master equation and the Redfield one at this stage. 

On the other hand, the Redfield equation predicts an evolution of the populations ruled by the following equation 
obtained by neglecting all the contributions coming from the coherences in Eq. (|51|) : 

f, ss (t) = 2^A 2 ]T|( S |5| S ')| 2 



x [F(e,e + E s ,-Eg) n(e + E s , - E s ) Vs , s ,(t) 
-F(e,e + Eg-E s ,)n(e + Eg-Eg,) Vss (t) } . 



(53) 
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This equation is the same as the master equation for the populations derived by Cohen- Tannoudji and co-workers in 
Ref. |f. 

We notice that important differences now exists between the population equation l|53(l obtained from the Redfield 
equation and the other population equation H42[) obtained from our master equation. Both equations describe the 
evolution of the populations as a random walk process in the spectrum. However, these processes are significantly 
different for Eqs. Ij42(l and Ij53(l . Let us focus on the evolution of the probability to be on a system state corresponding 
to the system energy E s . In both equations we see that for the loss contributions to the evolution coming from the 
jumps from an energy E s to an energy E s > , the density of states of the environment is modified by the energy E s — E s > . 
This is consistent with the Fermi golden rule applied to the total system and, thus, keeps the total energy constant. 
We care now on the gain contributions to the evolution. In our equation H42|) . we see that for these contributions 
due to jumps from the system energy E s i to E s , the density of states of the environment is modified by the energy 
E s i — E s . This is also consistent with the Fermi golden rule applied to the total system and, thus, keeps the total 
energy constant. However, for the Redfield equation we see that for the jumps from an energy E s i to an energy E s , 
the density of states of the environment is not modified by the energy E s > — E s , which is not consistent with the Fermi 
golden rule applied to the total system and does not keep the total energy constant. 

One can represent, for the Markovian case, the transitions described by Eq. (|42|) and i|53|) in a plane of the 
system energy versus the environment energy. In Figs. ^ and [2] we have depicted the energy exchanges described, 
respectively, by the Redfield equation and our equation in the Markovian limit for two different systems. Transitions 
between the system and the environment have to preserve the energy of the total system according to the Fermi 
golden rule and have therefore to occur along diagonal lines of the plane. One can see in Figs. ^ and [21 that only 
our equation satisfies this condition. The Redfield equation describes transitions that occur along a vertical line at 
constant environment energy and is therefore not consistent with energy conservation in the total system. This is 
acceptable if the environment is sufficiently large and has an arbitarily large energy. However, this is inadequate if 
the total energy of the system and the environment is finite as in nanoscopic systems, in which case our master should 
replace the Redfield equation. 

We can summarize as follows the differences between our quantum master equation (|38ll and the Redfield equation 
1)510. The derivation of both equations is based on the perturbative expansion of the total density matrix, but a 
specific form is imposed in each equation to the total density matrix [see Eqs. (|10[l and l|44|) ]. The consequence of 
this choice can be seen on the reduced density matrix of the system. In the Redfield theory, we have 

/5 S (t) = 5>)^(i) < S '|, ( 54 ) 

s,s' 

while, for our master equation, using Eq. I|12|) we have 

h(t) = Yl fde\ S )P ss ,(e;t){s'\. (55) 

s,s' J 

The system density matrix is related to the distribution functions according to 

(s\p S (t)\s') =w(<) = J deP ss ,{e-t) . (56) 

We see that, in our master equation, the matrix elements of the system density matrix are decomposed on the energy 
of the environment. This is not the case for the Redfield equation. The decomposition allows us to correlate the states 
of the system with the states of the environment. This is the main point of our master equation. The density matrix 
adopted for the Redfield equation cannot describe such correlations. In the Redfield equation, during the evolution, 
the environment is always in the same state while, in our master equation, the state of the environment is determined 
by the state of the system. As a consequence, we obtain a description which is consistent with energy conservation 
thanks to our master equation. 



III. APPLICATION TO THE SPIN-ENVIRONMENT MODEL 



In this section we consider a specific class of two-level systems interacting with an environment. The two-level 
system may be supposed to be a spin. An example is the spin-boson model in which the environment is a set of 
harmonic oscillators behaving as phonons [19| . 

The Hamiltonian of the spin-environment model we consider here is the following: 

fftot = j°z +Hb + \o x B . (57) 
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The eigenvalue equation of the system is 



Hs\s) = —v z \s) = sy|s) 



(58) 



where s — ±1. Like in Sec. ^ we first derive our master equation and then the Redfield equation in order to compare 
both equations. 



A. Using our master equation 

Let us now apply our master equation to the spin-environment model. 

In our theory and for a two-level system, the total density matrix becomes 

P(t) = ' 



n(H B ) 



P++(H B ;t)\- 



P+-(H B ;t) 



+P- + (A B ;t)\-)(+\+P—(H B ;t)\-)(-\ 

For the spin-environment model, our non-Markovian master equation l|36[l is given by 



P ++ (e;t) = -\ 2 P ++ (e;t) / de'F(e,e')n(e') / dr 



+X 2 n(e) I de'F(e,e')P—(e';t) J dr 



P—(e;t) = -A 2 P__(e;i) / de'F(e, e')n(e') / dr 



e i{e-e'+A)r + e -i(e-e'+A) 7 
e i(e-e'+A)r + & -i(e~e' +A)i 



e i(e-e'-A)T , e -i(e-e'-A)T 



+\ 2 n{e) / de' J F 1 (e,e')P++(e';t) J dr 



i(e— e'-A)i 



-i(e-e'— A)i 



P+_(e;t) = -iAP + _(e;i) 



\ 2 P+-{e-t) / de'F{e, e')n(e') / dr 



+A 2 n(e) / de'F(e,e')P- + (e';t) / dr 



e i(e-e'+A)T + e -i(e-e'-A)r 
e i(e-e'_A)r + e -i(e-e'+A)r 



and a further equation for P |^(e; t) given by the complex conjugate of Eq. (|(j2|) . 



(59) 



(60) 



(61) 



(62) 



We observe that the diagonal and off-diagonal elements of P ss ' (e; t) obey decoupled equations in the case of the 
spin-environment model. Therefore, the time evolution of the populations is independent of the time evolution of the 
quantum coherences. We now perform the Markovian approximation that consists of putting the upper bound of the 
time integral to infinity. Using Eq. (|37f) . we find 



P++(e; t) = 27rA 2 F(e, e + A) [n(e)P__(e + A; t) - n(e + A)P ++ (e; t)} 



(63) 



P—(e;t) = 2nX 2 F(e, e - A) [n(e)P ++ (e - A; t) - n(e - A)P__(e;i)] 
P+-(e;t) = -iAP + _(e;t) 



+iA 2 y de'P(e,e')P- 



2A 

! , - £ 2 - A 2 Ke )P+ - (e; t] + n{e)P - + (£ ; <)] 
-ttA 2 [n(e + A)F(e, e + A) + n(e - A)F(e, e - A)] P + _(e; t) 
+tt\ 2 n(e) [F(e, e + A)P_+(e + A; t) + F(e, e - A)P_+(e - A; *)] . 



(64) 



(65) 
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We notice that during the time evolution of the populations, the following quantity remains a constant of motion: 

P(e;t) = P ++ (e;t)+P__(e + A;t) =P(e;0) . (66) 
Accordingly, the difference of the populations defined as 

Z(e;i) = P ++ (e;t)-P__(e + A;i) (67) 

obeys the differential equation 

Z(e;t) = 2ir\ 2 [n(e)-n(e + A)}F(e,e + A)P(e;0) 



-2ttA 2 [n(e) + n(e + A)] P(e, e + A)Z(e; t) , 

the solution of which is given by 

Z(e; t) = Z(e; oo) + [Z(e; 0) - Z(e; oo)] e ~ 7Pauut 
with the asymptotic equilibrium value 

n(e) - n(e + A) 



Z(e; oo) 



n(e) + n(e + A) 



P(e;0) 



and the relaxation rate 



7Pauii = 2ttA 2 [n(e) + n(e + A)] P(e, e + A) 



(68) 
(69) 

(70) 

(71) 



Therefore, the populations relax to their asymptotic equilibrium values for each pair of energies e and e + A of the 
environment, keeping constant the initial distribution of the quantity P(e; 0). 

The time evolution of the distribution functions P±^(e;t) of the quantum coherences is more complicated because 
there is now a coupling between a continuum of values of the environment energy instead of only two values. Accord- 
ingly, the distributions of quantum coherence is ruled by a couple of two integro-differential equations, instead of an 
ordinary differential equation. 

B. Using the Redfield equation 

For the spin-environment model, the non-Markovian and Markovian Redfield equations can be derived from Eqs. 
gEJ and J5TJ. 

Using Eq. I|48(l , the non-Markovian Redfield equations here write 



f)++(t) = -X 2 v++(t) / de'n(e')F(e,e') / dr 



+A 2 77 (t) / de'n(e')F(e, e') J dr 



»)__(t) = -A 2 77 (t) / de'n(e')F(e, e') / dr 



-\ 2 V++ (t) I de'n(e')F(e,e') J dr 



i(e-e'+A)T , & — i(e-e'+A)7 



e i(e-e'-A)r _|_ g-i(e-e'-A)r 



e i(e-e'-A)r _|_ g -i(e-e'-A)T 



,t(e-e'+A)r _|_ g -t(e-e'+A)7 



(72) 



(73) 



»)+-(*) = -»A>)+-W 



A 2 »7 + _(t) / de'n(e')P(e, e') / dr 



+A 2 r/_+(t) I de'n(e')F(e,e r ) J dr 



J(e-e'+A)T + e _i( 6 _ e '-A)T 



e i(e-e'-A)r + e -i(e-e'+A)r 



(74) 
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and a further equation for rj |_(t) given by the complex conjugate of Eq. (|74|l . 

Here again, there is a decoupling between the time evolutions of the populations and of the quantum coherences. 

Taking the Markovian approximation by replacing J* into J °° and using Eq. (|37|l . we get the Markovian Redfield 
equations for the spin-environment model: 

T) ++ (t) = 27rA 2 [n(e-A)F(e,e-A)r;__(t) 

-n(e + A)F(e,e + A)r )++ (t)] , (75) 

fi—(t) = 2ir\ 2 [n(e + A) F(e, e + A) 77++ (i) 

-n(e-A)F(e,e-A)r ? __(t)] , (76) 



i7+_(t) = -<A»7 + _(t) 

2 

-ttA 2 [n(e + A)F(e, e + A) + n(e - A)F(e, e - A) 



of,, , 2A 
+J A 2 y de'n(e')F( e , £ ')^^37)2-rA^ (*)+ »?-+(*)] 



x [v+-(t)-V-+(t)] ■ 

(77) 

The populations of the two-level system are controlled by the z component of the spin defined as the difference 

ZRedfleld(t) = V++(t) - V — (t) . (78) 

According to the Markovian Redfield equations (|75|l and l|76|) , the z component obeys the differential equation 

iRcdfioid = 2 7 rA 2 [n(e-A) J F 1 (e,e-A)-n(e + A)F(e,e + A)] 

-2ttA 2 [n(e - A)F(e, e - A) + n(e + A)F(e, e + A)] z Rod ficid ■ 

(79) 

Its solution is given by 

ZRedfieldW = 2 Rcd field(oo) + [z Rc dficld(0) - 2Rcdficld(oo)] e~ 1Redtieldt , (80) 

with the equilibrium value 

n(e - A)F(e, e - A) - n(e + A)F(e, e + A) 
ZRcdficld(00) - n( e -A)f( £!e -A) + n( e + A)F( e , e + A) ' (81) 

and the relaxation rate 

7Rcdficid = 2^A 2 [n(e - A)F(e, e - A) + n(e + A)F(e, e + A)] . (82) 

We notice that the rate predicted by the Redfield equation coincides with the one predicted by our master equation 
only in the limit A — 0. A more important difference appears in the asymptotic equilibrium values for the z component 
of the spin predicted by both equations. These differences find their origin in the problem of conservation of energy 
with the Redfield equation, as explained above. Comparison with numerical data will confirm this explanation in a 
following section in the case of the spin-GORM model. 

IV. APPLICATION TO THE SPIN-GORM MODEL 

In order to confront our master equation and the Redfield equation with numerical data and test their respective 
domains of validity, we now apply our theory to a specific class of two- level systems interacting with an environment, 
for which the environment operators are Gaussian orthogonal random matrices (GORM). We call this model the 
spin-GORM model and its detailed properties will be described elsewhere |20|. 
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The system is a two-level system, while the environment is supposed to be a system with a very complex dynamics. 
Here, the term complex is used in a generic way. The complexity can come, for example, from the fact that the 
corresponding classical system is chaoticlike in a quantum billiard or for the hydrogen atom in a strong magnetic 
field |2lL l22j . It can also come from a large number of coupling between states in an interacting many-body system 
like those appearing in nuclear physics |22j or in systems of interacting fermions like quantum computers |22j. A 
well-known method, developed by Wigner in the 1950s, for modeling the energy spectrum of a complex quantum 
system containing many states interacting with each other, consists of assuming that their Hamiltonian is a random 
matrix |2.'l l24l |25| . Here, we suppose that the Hamiltonian of the environment is a Gaussian orthogonal random 
matrix (GORM). The interaction between the spin and the environment is given by a coupling operator which is the 
product of a system and a environment operators. The latter is also represented by a GORM because of its complex 
interaction with the many degrees of freedom of the environment. Such random-matrix models have recently turned 
out to be of great relevance for the discussion of relaxation and dissipation in quantum systems 0, |2|| [23, l28l |2^| . 

The spin-GORM model can therefore be considered as a particular case of the spin-environment model in which 
the environment operators are GORM. The Hamiltonian of the spin-GORM model is thus given by 

#tot = +H B + Xa x B (83) 



where the Hamiltonian of the environment is 



and the environment coupling operator by 



2 



Hb = -= X , (84) 



B = -== X' . (85) 



X and X' are two statistically independent ^ x ^ random matrices belonging to the Gaussian orthogonal ensemble 
(GOE) of probability density 

p{X) =C exp (-a^TrX 2 ) , (86) 

with a-x = \ and a normalization constant C. N/2 is the number of states of the environment. The off-diagonal 
and diagonal elements of X are independent Gaussian random numbers with mean zero and standard deviations 
CToff-diag = 1 an d CTdiag = V2 respectively. 

In the limit N — > oo, the density of states of the environment gets smooth and can be calculated by an average 
over the random-matrix ensemble 



N/2 

n(e) = ^2S(e-E b ). (87) 



6=1 

It is known that the GOE level density is given by the Wigner semicircle law: 



1 ,2 it 1,1 ^ 1 



n(e) = <j —VI c if | e |< 2 (88) 
t0 if |e|>|. 

The random matrices are normalized so that the level density of the environment has a width equal to unity. To 
simplify the notations in the following, we use the convention 



For the following, we also need to evaluate the function F(e, e') for our random-matrix model. For this purpose, we 
need the random-matrix average of the quantity l|28|). Since, in the GOE, we have that 
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we find that 



Tr B S(e- H B )BS(e' - H B )B 
Y l S(e-E b )6(e>-E b ,)\(b\B\b>)\< 

b.b' 



^5(e - E b )6(e' - E b )\(b\B\b)\* + £ S(e - E b )5(e' - E b ,)\(b\B\b'] 



b b^b 



= i^E ^ - E ^' E ») + ^ E 6 ( e - - E » ) 

b b^V 

±= S (e ~ E b )S(e' -E b ) + ^-J2 ^ e " Eb ^ € ' ~ Eb ' ) 



8N^ v ' v ' 8N< 

b 6,6' 

~ ^( e - e ')n(e) + ^n(e)n(e'). (92) 

In the limit TV — > oo, the first term becomes negligible in front of the second term so that the comparison with Eq. 
(J2HJ shows that for the spin-GORM model, 

i?( ee ')~J_ (93) 
v ' ; 8N y ' 

The total system contains N states. The unperturbed density of states of the total system is schematically depicted 
in Fig. |3 for A = 0. The model has different regimes whether the splitting A between the two levels of the spin is 
larger or smaller than the width of the environment level density. The spin-GORM model can describe a large variety 
of physical situations. In the present paper we focus on the perturbative regimes (A -C H B ). When A is larger than 
the width of the semicircular density of states of the environment, we are in a highly non-Markovian regime. The 
dynamics of the system is faster than that of the environment. On the other hand, when A is smaller than unity, we 
are in a Markovian regime because the dynamics of the environment is much faster then the one of the system. 

Now, we apply our master equation and the Redfield equation to the spin-GORM model in both their Markovian 
and non-Markovian versions. 



A. Using our master equation 

We now apply our master equation (|36[) to the spin-GORM model. Using Eqs. (|60|l - l|t)2|) . (JSSJ, and (|93|) . we get 

the non-Markovian equations 



P++(e;t) / de'M- r 



A* 

a 2 r\ 



7T V 4 



A 



e 2 / cfe'P__(e';i) 



P__(e;t) = P— (e;i) / de\ - - e 



A 2 1 



7T V 4 



/ £ 2 / de ' P++{e ' ;t) 



sin(e — e' 


+ A)t 


(e-e'- 


fA) 


sin(e — e' 


+ A)t 


(e-e'- 


f A) 


sin(e — e' 


- A)t 


(e-e'- 


-A) 


sin(e — e' 


- A)i 


(e-e>- 


-A) 



(94) 



(95) 



P+-(e;t) = -iAP+_(e;t) 
A 2 

A 2 / 1 



^P+-(e;t) I de\/~-e' 2 / dr 



2?r V 4 



e 2 / de'P_ + (e';t) / dr 



J(e-e'+A)r _j_ e -i(e-e'-A)7 



e i( £ -e'-A)r + e -i(e-e'+A)T 



(96) 
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We notice that the equations for the populations are decoupled from the ones for the quantum coherences. 



Performing the Markovian approximation and using lim 4 



= ttS(uj) and Eq. (|37Jl . we get 



/'++(<:')= -X\j\ - (e+ A) 2 P++ {e- 1) + AV J - e 2 P__(e + A; t) , 



(97) 



P_ !<:/): -A 2 \/i-( e -A) 2 P__( e ;i) + Ayi- e 2 P++ ( e _A ; t) 



(98) 



P+_(e;t) = -iAP+_(e;i) 
. A 2 



l-e' 2 P + -{e-t) + J\-^P_ + (e'-t) 



A 2 
2 



J-( e + A) 2 + ^-( e -A) 2 



A 2 1 



— W - - e 2 [P_ + (6 + A; t) + P_ + (e - A; t)} 



(99) 



where the expressions and integrals over energy extend over the interval of definition of the level density n(E) and of 
the distributions P ss /(E;t) which is always —1/2 < E < +1/2, E being the argument of these functions. 
We now focus our attention on the evolution of the populations. 

We see from Eqs. (19 7|) and (|98|l that the transitions conserve the total energy of the system and environment so 
that the transitions occur between the only two energies e and e + A of the environment. As a consequence, the 
quantity 



P(e; t) ee P ++ (e; t) + P__(e + A; t) = P(e; 0) 



(100) 



is a constant of the motion for each energy e of the environment, as already noticed with Eq. I|66|) . The difference 
(I67|l of populations obeys the differential equation (|68() . If the initial distributions P(e'; 0) and Z(e'; 0) are Dirac delta 
distributions centered on the initial energy e: 



P(e';0) = 8(e'-e), 
Z(e';0) = J(e'-e) z Pauli (0) 

The z component of the spin defined as 

zpauli(*) = J de Z(e;t) 
obeys the same differential equation as the distribution Z(e; t), 

ipauli(i) = A 2 



(101) 
(102) 



(103) 



1 2 

4~ £ 



i -(^ + A) 2 



-A 



1 



1 



(e + A) 2 



ZPauli(t) 



The solution of Eq. (|104|) is given by 

ZPauli(t) = ^Pauli(oo) + [z Pauli (0) - Zp a uli(oo)] e~ 7Pauli * , 

with the asymptotic equilibrium value 



^Pauli(oo) = 



i - f 2 
4 fc 



I-(^ + A) 2 



(104) 



(105) 



(106) 
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and the relaxation rate 



7Pauli 



A 2 



(107) 



With the convention Q89p. the expressions are nonvanishing only over the interval of definition of their argument. 
Figure 0] helps us to represent the different values that take Eqs. (|106fl and (|107fl in the space of the environment 
energy e and of the splitting energy A of the two-level system. 



B. Using the Redfield equation 

For the spin-GORM model, the correlation function (|46(l can be calculated by performing a GOE average. Using 
the level density H88|) and the value l|93l) . the microcanonical correlation function l|5U|) becomes 

a(r,e) = / de' n(e') F(e, e') e l{e ~ e ' )T 



it V 4 87V 



At 



e ltT , (108) 



in the limit N — > oo, where J\{t) is the Bessel function of the first kind. 

Therefore, using the Redfield equation (|72() - (|74|l and the microcanonical correlation function of the spin-GORM 
model, we get 



»?++(*) = -A 2 ?7++(0 / drcos[(e + A)r 
Jo 



2r 
•Mf) 



+A 2 / drcos[(e- A)r]4H - ( 109 ) 

Jo 2r 



77 (t) = —A 2 77 (t) / drcos[(e - A)r 

Jo 



Mi) 

It 



+.\- //_+(/)/ drcos[(e + A)r]^ , (110) 



!)+_(«) - -iA?7 + _(t) 

-A 2 7?+-(t) / dre jAT cos(er) 



- .. ,/, / ,/-, 'A: ^- Mj) 

2r 



+ A 2 ,7. .(/) f dTe~ lAT cos(eT)^l . (Ill) 



These are the non-Markovian Redfield equations for the spin-GORM model. 

The Markovian Redfield equations for the spin-GORM model take the following forms: 



'/++!/ 1 = -A 2 ^-(e + A)^ ++ (t) 



+\\/\-(e-£) 2 V—V)> (112) 



'/ (') = -A 2 ^/i-(6-A)^__(t) 



+ A 2 A /i-( e + A)^ ++ (i), (113) 
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-iA v+ -(t) 
A 2 



' < ' - A 2 
x [rn._(t) + »?_+(<)] 



2 



i -(^ + A)2 



(114) 



We focus on the evolution of the populations. 

The population of the two- level system is controlled by the z component of the spin by Eq. Ij78[l . According to the 
Markovian Redfield equations (|112J) and (|113|l . the time evolution of the z component is given by Eq. (|80|1 with with 
the equilibrium value 



ZRedfield(oo) = 



i-(e-A)2- ,A-(e + A) 



i-(e-A)2+ l-( e + Ay 



(115) 



and the relaxation rate 



\2 

7Rcdficld — A 



\-{e-AY + \J\-{e + Ay 



(116) 



Figure 0] depicts the different regimes predicted by this equation in the space of the environment energy e and the 
splitting energy A. 

We observe that in the limit A — > 0, where the energy scale of the system is much smaller than the one of the 
environment, both our master equation and the Redfield equation predict a similar value for the relaxation rate. 
However, differences appear for the value of the asymptotic value of the z component of the spin. As we explained 
here above, the reason is that our master equation is consistent with energy conservation in the total system albeit 
the Redfield equation is not. This problem spoils the applicability of the Redfield master equation if the environment 
is not arbitrarily large as we shall see in the next section. 



V. NUMERICAL RESULTS AND DISCUSSION 



The purpose of the present section is to compare the different master equations with exact numerical results 
obtained for the relaxation of the z component of the spin due to the interaction with its environment in the spin- 
GORM model. The initial condition of the spin is always the state |+). The environment is always taken in a 
microcanonical distribution at a given energy e. The width of the energy shell of this microcanonical distribution is 
always equal to Se = 0.05. 

A general comment is here in order concerning the applicability of a master equation to a quantum system with 
a discrete energy spectrum. Indeed, beyond a time longer than the Heisenberg time (which is defined as the level 
density of the total system) , quantum beats and recurrences appear due to the discreteness of the energy spectrum. 
Only, the decay before the Heisenberg time can be compared with the prediction of a quantum master equation. It 
turns out that the further condition AA 2 > 1 should also be satisfied, which requires that the coupling parameter 
should not be too small with respect to the mean level spacing which goes as 1/A. If this condition is not satisfied 
(AA 2 < 1) the time evolution of individual systems present large quantum oscillations which widely deviate from the 
prediction of the master equation. On the other hand, if A A 2 > 1, the deviations with respect to the predictions of 
the master equation are smaller than the signal itself and tend to decrease as A — > oo [2(j. In the limit A — > oo, the 
decay curve which is the solution of the master equation is approximately followed by a majority of realizations of the 
process by individual systems. In the figures given here, these deviations are not seen because of an averaging of the 
signal over x — 10 individual systems. Besides the condition A A 2 > 1, the coupling parameter should also be small 
enough to justify the perturbative treatment, typically A < 0.3. 

The equations we are comparing are the following: 

• The von Neumann equation describes the z component of the spin using Eq. (J2J) with the Hamiltonian i|83[) . 
Averaging is carried out with \ realizations of the GORM Hamiltonian. This calculation does not involve any 
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approximation and, therefore, gives the exact solution of the problem. All the following equations will be compared 
to this one. 

• The most general non-Markovian version of our master equation (j^fijl using Eqs. 1)94(1 and (|95[1 . which we refer to 
as the Pauli non-Markovian (NM) equation. 

• The Markovian version l|38|) of our master equation, which we refer to as the Pauli Markovian (M) equation. For 
the spin-GORM model, this equation is given by Eqs. l|9Tjl and JHEl and its solutions by Eqs. H105|) - (|107jl . 

• The Redfield non-Markovian equation (gSJ is given by Eqs. (|ll)9f> and (|110fl for the spin-GORM model. 

• The standard Redfield Markovian equation l|51|l is given by Eqs. (|112f> and i|113|) and its solutions by Eqs. I|8U|) 
with Eqs. (|TT5j) and (jTToj) . 

The results of the numerical calculation of the time evolution of the z component of the spin are depicted in Figs. 
l5ll(Jl for different regimes of the spin-GORM model, i.e., for different values of the energy splitting A of the two-level 
system as well as of the environment energy e. In all the cases, the coupling parameter is equal to A = 0.1. 

Figures 03 and |H1 depict the global relaxation of the z component of the spin for increasing values of the energy 
splitting A. In accordance with what we argued before on theoretical grounds, we see in these figures that the larger 
the energy splitting A of the system is, the bigger is the difference between the Redfield and our master equation. We 
also see that our equation always fits very well with the exact von Neumann equation, which is not the case of the 
Redfield equation. As argued before, this is due to the fact that the Redfield equation does not take into account the 
changes in the energy distribution of the environment induced by the system transitions. When the system energy is 
very small this makes almost no difference, but when it increases, this has to be taken into account and our master 
equation becomes necessary. 

In particular, a large discrepancy happens for the solution of the Redfield equation in Fig. [S] although the solution 
of our master equation continues to be in agreement with the exact time evolution. This can be understand with Eqs. 
(|l(J6fl and (|115fl for the asymptotic equilibrium values of the z component, which predict, respectively, 

^Pauli = 0, (117) 
^Redfield = — 1 ; (H8) 

for A = 0.5 and e = —0.25. The discrepancy of the Redfield equation finds its origin in the violation of energy 
conservation between the system and its environment by this equation. The Pauli equation has the advantage of 
allowing a correct energy exchange between the spin and its environment, which is crucial for obtaining the correct 
asymptotic equilibrium value of the z component. 

We can also notice in Fig. [3] and Figs. 17181 that, on the global time scale, the non-Markovian and Markovian 
equations are very close to each other. But if we look on a shorter time scale, we see in Figs. [H] and [!|] small 
differences between the non-Markovian and Markovian equations in the early stage of the decay. The solutions of 
the non-Markovian equations are in best agreement with the exact time evolution and present a non exponential 
early decay on the time scale of the environment correlation time (r corr ~ 10). In contrast, the solution of the 
Markovian equation immediately enters in an exponential decay and, thus slightly deviates from the exact solution. 
This observation concerns both the Redfield and Pauli equations. This suggests that, as explained in Refs. [TTL IT^ . 
a slippage of initial conditions is required for both Markovian equations in order to avoid this small early-decay 
discrepancy. We also observe that this discrepancy decreases with the energy splitting A. This is expected since the 
Markovian approximation is valid if the time scale of the system dynamics ^ is longer than the environment time 
scale (r corr ~ 10). Therefore, the smaller A is, the better is the Markovian approximation. If one wants to make a 
correct description of the system dynamics on a time scale of order r corr , the non-Markovian equations should be used 
(or the Markovian equations should be supplemented by a slippage of initial conditions 0,0). We also see, in Fig. 
El that the Pauli non-Markovian equation gives better results than the Redfield non-Markovian equation not only on 
long time scales but also on short ones, even if their solutions essentially coincide on a very short time of order r corr 
but not more. Again, when the energy of the system is to small to affect the density of state of the environment, this 
difference between our non-Markovian equation and the Redfield non-Markovian equation disapears (see and 0). 

We see in Fig. ^| the relaxation of the z component of the spin in a highly non-Markovian regime. The energy 
difference between the two levels is here much larger than the width of the level density of the environment. The 
Markovian equations are not plotted here because they describe a constant value equal to unity for all times and, 
therefore, completely miss the dynamics. The whole spin dynamics happens on a time scale of order r corr . We also 
see that, in this highly non-Markovian regime, there is almost no difference between the non-Markovian Redfield 
equation and our master equation on the short time scale that we plotted. The non-Markovian Redfield equation as 
well as our non-Markovian master equation continue to fit with the exact dynamics even on longer time scales. The 
special structure seen in Fig. ^]can be understand by using Eqs. (|109fl and l|110fl . Indeed, the curve is the result of 
some time integrations of the Bessel function Ji(t/2) divided by t and modulated by cosAi. Since the modulations 
of cos At have a period ^ shorter than the decay time r corr ~ 10 of the Bessel function Ji(i/2), a shape reminiscent 
of a Bessel function only appears as an envelope of the oscillations of the decay curve. 
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We conclude that, as expected from theoretical arguments, our master equation gives excellent predictions, especially 
in situations where the system energy is greater than or of the same order of magnitude as the typical energy scale of 
variation of the density of states of the environment. We also notice that for non-Markovian dynamics that happen 
on a time scale of order r corr , the non-Markovian Redfield equation give the same result as our master equation for 
short time scales. But for longer time scales our equation is the only one that correctly describes the dynamics. 

VI. CONCLUSIONS 

We derived in this paper a master equation to study the dynamics of a quantum system interacting with its 
environment. This equation is obtained by a perturbative expansion with respect to the coupling parameter between 
the system and its environment. Our equation is more general than the previously obtained perturbative master 
equations because our equation explicitly takes into account the exchange of energy between the system and its 
environment. This effect is important when the density of state of the environment varies in a significant way on 
energy scales of the order of the system energy scales. 

We showed how the well-known master equations of the literature can be derived from our equation (|36|l by 
performing different types of approximations. Our equation reduces to the non-Markovian Redfield equation (|48ll if 
one neglect the changes in the density of states of the environment induced by the system. Moreover, we showed 
that by performing the Markovian approximation on our equation and neglecting the coherence contribution to the 
population dynamics, we get a Pauli-type equation (|42|l for the total system (system + environment) that describes 
the time evolution in terms of distributions defined on the energy spectrum of the environment. When one neglects the 
changes in the energy distribution of the environment, the Markovian version of our equation (that is now equivalent 
to the Redfield equation) reduces to the master equation derived by Cohen- Tannoudji and corworkers in Ref. Q- 

We have applied our equation to a two-level system interacting with a general environment (spin-environment 
model), especially, in the case where the environment operators are random matrices (spin-GORM model). In this 
case, we have carried out numerical simulations of the spin-GORM model (for which one can computed the exact 
solutions) that show the greater accuracy of our master equation with respect to the other well-known master equations 
in Markovian and non-Markovian situations. 
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FIG. 1: Schematic representation of the energy exchanges described, respectively, by the Redfield and by our master equation 
in the Markovian limit for a two-level system model, in the plane of the system energy Es versus the environment energy Eb- 
The energy splitting between the two levels of the system is denoted by A. The energy spectrum of the system is discrete 
(two levels) while the one of the environment is a quasicontinuum represented by the density of states given by the Wigner 
semicircle law 188H of width equal to unity. The total energy of the system is given by E — Es + Eb, which corresponds to the 
diagonal line. The initial condition is denoted by two empty superposed circles. We see that transitions preserving the total 
energy have to occur along the diagonal line E — Es + Eb- Doing this, they satisfy the Fermi golden rule for the total system. 
One can see that only our master equation satisfies this condition (dotted transition lines). The Redfield equation describes 
transitions that occur along a vertical line at constant environment energy and is therefore wrong (dashed transition lines). 
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FIG. 2: Generalization of the previous Fig. 0for the case where the system has more than two levels (here four levels). One 
can see that the system levels (horizontal lines) that do not intersect the total energy diagonal line E = Es + Eb within the 
environment energy spectrum delimited by the sparse-dotted vertical lines do not participate in the dynamics. 
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FIG. 3: Schematic representation of the smooth density of states (DOS) of the unperturbed spin-GORM model (A = 0) for 
different values of the energy splitting A of the spin. The horizontal axis is the environment energy e while the vertical axis is 
the energy splitting A. In the lower and central parts, the splitting A is smaller than the width of the environment DOS. In 
the upper part, the splitting A is larger than the width of the DOS. 
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FIG. 4: Schematic representation of the different regimes of the spin-GORM model for situations where the initial state of the 
spin is s = +1, in the plane of the environment energy e versus the spin energy splitting A. The different regions correspond 
to different values of the functions n(e), n(e + A), and n(e — A), where n(E) denotes the DOS defined by the semicircle law 
11881 . One can take a value of the environment energy anywhere between e = — | and e = |, where n(e) 7^ 0. In region 1: 
n(e + A) / and n(e - A) = 0. In region 2: n(e + A) = and n(e - A) / 0. In region 3: n(e + A) / and n(e - A) ± 0. In 
region 4: n(e + A) = and n(e — A) = 0. This figure is useful to evaluate equations such as Eqs. I|l()fi|l and l|l()7fl and l|115|l 
and 11161 . 
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FIG. 5: Relaxation of the z component of the spin for the spin-GORM model for a very small spin energy splitting A = 0.01 
with A = 0.1, e = 0.25, N — 2000, and \ = 10- The exact curve is given by integrating the von Neumann equation, which is 
compared with the solutions of the Pauli and Redfield Markovian (M) and non-Markovian (NM) equations. We see that all 
the perturbative equations give similar results in the present case. 




FIG. 6: Relaxation of the z component of the spin for the spin-GORM model in the same conditions A = 0.01, A = 0.1, 
e = 0.25, N = 2000, and \ = 10 as in Fig. |S]in order to show that, on a short time scale of the order of the correlation time 
of the environment t CO tt — 10, the non-Markovian equations (denoted by NM) describe very accurately the dynamics although 
the Markovian ones (denoted by M) is exponential and deviate from the exact behavior. On a longer time scale (much longer 
than Tcorr), the solutions of the Markovian equations join those of the non-Markovian equations. 
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FIG. 7: Relaxation of the z component of the spin for the spin-GORM model for a small spin energy splitting A = 0.1 with 
A = 0.1, e = 0.25, N — 2000, and \ = 10. The exact solution of the von Neumann equation is compared with the solutions of 
the Pauli and Redfield Markovian (M) and non-Markovian (NM) equations. We see that our master equation (Pauli NM) gives 
the best results and that the solutions of the Markovian equations remain very close to those of the non-Markovian equations. 
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FIG. 8: Relaxation of the z component of the spin for the spin-GORM model for an intermediate spin energy splitting A = 0.5 
with A = 0.1, e = —0.25, N = 2000, and \ — 10. The exact solution of the von Neumann equation is compared with the 
solutions of the Pauli and Redfield Markovian (M) and non-Markovian (NM) equations. Here, we see that the Redfield equations 
give completely wrong results after a short time. The Pauli equations give much better results than the Redfield ones. 




FIG. 9: Relaxation of the z component of the spin for the spin-GORM model in the same conditions A = 0.5, A = 0.1, 
e = —0.25, N = 2000, and x — 10 as in Fig. |H] We focus here on the short time dynamics in order to see that only the 
non-Markovian equations (NM) reproduce the initial behavior of the system which is not the case for the Markovian equations 
(M). After r corr ~ 10, the non-Markovian Redfield equation (Redfield NM) becomes wrong but our master equation (Pauli 
NM) is still valid. On a longer time scale, the solution of the Markovian version of our master equation (Pauli M) joins the one 
of the non-Markovian version of our equation (Pauli NM) . 
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FIG. 10: Relaxation ol the z component of the spin for the spin-GORM model for a large spin energy splitting A = 5 with 
A = 0.1, e = 0.25, N = 2000, and \ = 10. The exact solution of the von Neumann equation is compared with the solutions of 
the Pauli and Redfield non-Markovian (NM) equations. We are in a highly non-Markovian regime. The Markovian equations 
are not plotted here because their solutions are a constant equal to unity at all times and therefore miss the whole dynamics. 
We see that in this highly non-Markovian regime and on the short time scale r CO rr there is almost no difference between the 
non-Markovian Redfield equation (Redfield NM) and our master equation (Pauli NM). 



